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AS DISCRETE PAINLEVE EQUATIONS. 



Alphonse P. Magnus 

Abstract. We consider orthogonal polynomials p n with respect to an ex- 
ponential weight function w(x) = exp(— P(x)). The related equations for 
the recurrence coefficients have been explored by many people, starting 
essentially with Laguerre [49], in order to study special continued frac- 
tions, recurrence relations, and various asymptotic expansions (G. Freud's 
contribution [28,56]). 

Most striking example is n = 2tw n + w n (w n+ i + w n + w n -i) for the 
recurrence coefficients p n +i = %Pn — w n p n -\ of the orthogonal polynomials 
related to the weight w(x) = exp(— 4(tx 2 + x 4 )) (notations of [26, pp. 34- 
36]). This example appears in practically all the references below. The 
connection with discrete Painleve equations is described here. 



1. Construction of orthogonal polynomials recurrence 

coefficients. 

Consider the set {p n }o° of orthonormal polynomials with respect to a 
weight function w on (a part of) the real line: 

Pn(x)p m (x)w(x) dx = 8 n , mi n,m = 0,l,... (1) 

The very useful recurrence formula has the form a\pi{x) = (x — b )p 0} 

a n+1 p n+1 (x) = (x - b n )p n (x) - a n p n _ 1 (x), n = l,2, ... (2) 

The connection between the nonnegative integrable function w and the real 
sequences {a n }^°, {& n }o° is of the widest interest. It is investigated in combina- 
torics [9,10,52], asymptotic analysis [11,12,28,30,53-55,61,77-79], numerical 
analysis [15,21,31,32,51,57], and, of course, Lax-Painleve-Toda theory (all 
the other references, excepting perhaps Chaucer). 

However, the connection may seem quite elementary and explicit: 

1. From w to a n and b n : let /i^ = x k w(x) dx, k = 0,1,... be the 
moments of w, then 

a n= (T r x 2 ,n = 1,2,. .. b n = — — ,n = 0,l,... (3) 

\tl n ) tin+1 tin 
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where H, n and H, n are the determinants 



n = 0,1,..., #„ = 1,H = 0. 
2. From the a„'s and 6„'s to w: consider the Jacobi matrix 



(4) 



a i then, remarking that (2) may be writ- 

a i & i a 2 ten xp n (x) = a n p n _i(x) + b n p n (x) + 

a 2 ^2 a 3 a n _|_ip n+ i (x), (Remark that the coefficients ap- 

pear now with positive signs, a very elementary 
fact, but always able to puzzle weak souls), or x p(x) = Jp(x), where p(x) 
is the column (infinite) vector of p (x) } pi(x) } . . . , so that the expansion 
of x k p n (x) in the basis {po,pi, . . . } is x k p n (x) = J^(«/ )n,tPt(%)- As we 
are dealing with a basis of orthonormal polynomials, this means that 
(J k ) Uj £ is the £ th Fourier coefficient of x k p n (x) 1 (where row and column 
indexes start at zero), i.e., 
polynomial p } one has 



(J )n,t = J^oo x k p n (x)pt(x)w(x) dx. For any 



(P(J)) 



nl 



p(x)p n (x)pi(x)w(x) dx. 



(5) 



In particular, («7 fc ) ,o = f^k/ Vo- 
Direct numerical use of (3) is almost always unsatisfactory, for stability 
reasons. Together with various ways to cope with the numerical stability 
problem [31,32], compact formulas or equations for the recurrence coefficients 
a n and b n in special cases have been sought. Of course, quite a number of exact 
solutions are known (appendix of [14]) and new ones are steadily discovered 
(Askey- Wilson polynomials and other instances of the new u q— calculus" [4, 
35,36,47]). The trend explored here is essentially related to weight functions 
satisfying simple differential equations. Main innovators were Laguerre in 
1885 [49], Shohat in 1939 [77], and Freud in 1976 [28]. Their contributions 
are now described in reverse order. On Geza Freud (1922-1979), see [29,68,69]. 



2. Definition of Freud's equations. 

In order to reduce the technical contents of what will follow, only even 
weights: V real x, w(—x) = w(x) will be used. Then, the orthogonal poly- 
nomials are even or odd according to their degrees; equivalently, all the b n } s 
in (2), in (3), and in the Jacobi matrix J vanish. Let a be the sequence of 
coefficients {a i5 a 2 , . . . }. 

Theorem 1. Let P be the real even polynomial P(x) = c x 2m + Cix 2m ~ 2 + 
' ' ' + c m; with c > ; and let ai,a 2 , ... be the recurrence coefficients (2) of 
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the orthonormal polynomials with respect to the weight w(x) = exp( — P(x)) 
on the whole real line. Then, the Freud's equations 

F n (a):=a n (P'(J))„,„-i=n, n = l,2,... (6) 

hold, where ( ) n>n -i means the element at row n and column n — 1 of the 
(infinite) matrix P'(J). A finite number of computations is involved in each 
of these elements. The row and column indexes start at 0. We could as 
well take the indexes (n — as J and any polynomial function of J are 

symmetric matrices. 
Remark. One has 

F 1 (a) = (JP'(J)) 0fi . (7) 

Indeed, the first element of the first row of JP'(J) is a\ times the element 
(1,0) of P'(J), i.e., F 1 (a). 
Examples. 



w(x) = exp(— x 2 ) F n (a) = 2a 2 n , (8) 
w(x) = exp(— ax 4 — fix 2 ) F n (a) = 4aa 2 l (a 2 l _ 1 + a 2 n + a 2 n+l ) + 2/3a^, (9) 
w(x) = exp(-x 6 ) F n (a) = §a 2 a (a 2 a _ 2 a 2 a _ x + a A a _ x + 2a 2 a _ x a 2 a + a A n 



,2 2 ,o22 i 4 i 2 2\ 

+ + 2a n a n+1 + a n+1 + a n+1 a n+2 ) } 



(10) 



The elementary Hermite polynomials case is of course immediately recovered 
in (8). (9) and (10) were used by Freud [28] in investigations on asymptotic 
behaviour, see [11,53,56,57,67,69] for more. The rich connection of (9) with 
discrete and continuous Painleve theory will be recalled later on. For a much 
more general setting, see [1, Theor. 4.1]. 

Proof of (6). We consider two different ways to write the integral 

/>oo 

In = / p , n (x)p n _i(x)exp(-P(x))dx. First, let p n -i(x) = -K n _ x x n ~ x + ... 



and p n (x) = ir n x n + . . . . From (2) (with n — \ instead of n), 7r n = 7r n _i/a n ; so, 

n 

p' n (x) = mr n x n ~ l + . . . = — Pn-i ( x ) J t~ a polynomial of degree ^ n — 2. As p n -\ 

a n 

n 

is orthogonal to all polynomials of degree ^ n — 2, what remains is I n = — , 

a n 

as the orthonormal polynomials have a unit square integral. Next, we per- 

/oo 
Pn(x) \pn-i (x) exp( — P(x))] dx = 
-oo 

P' (x)p n (x)p n -\ (x) exp( — P(x)) dx } using the orthogonality of p n and p' n _ 1 



of degree < n — 1 < n, and the latter integral is (P'(«7)) njn _i , according to 
the spectral representation (5). □ 
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3. Freud's equations software. 



The interested reader (if there is still one left) will find in the "soft- 
ware" part of http://www.math.ucl.ac.be/~magnus/ three FORTRAN pro- 
grams related to Freud's equations, freudl.f asks for an exponent m and 
puts data associated to the Freud's equations for the file 

freud00m.dat. freud2.f produces a readable (sort of) listing out of such a 
data file. Finally, freud3 .f reads the coefficients of a polynomial P and com- 
putes in a stable way a sequence of positive recurrence coefficients a i} a 2 , . . . 
associated to the weight exp( — P(x)) on the whole real line. Actually, weight 
functions \x\ a exp( — P(x)) } with P even, on the whole real line, and 
x 13 exp( — Q(x)) on the positive real line are considered simultaneously (P(x) = 



We started from the solution (3), and built the equation (6) afterwards, 
when the a n 's are recurrence coefficients of orthogonal polynomials associated 
to the weight exp( — P(x)). 

We try now to discuss the general solution of (6). With P even of degree 
2m, the n th equation of (6) involves a n _ m+ i, . . . , a n+m _i (see [56,57] or the 
examples (8) (9), (10)), so that the general solution depends on 2m — 2 
arbitrary constants, for instance a no _ m+ i, . . . , a no+m _ 2 . I shall only solve the 
case a_ m+ i = • • • = a = 0, and investigate how the solution depends on the 
m — 1 initial data a i} . . . , a m _i: 

Theorem 2. With P(x) = c x 2m + cix 2m ~ 2 + • • • + c m , the solution of (6) 
with a_ m+ i = • • • = a = is (3), with (4) ; where the fj, 3 's of (4) satisfy the 
linear recurrence relation of order 2m 

2mc fi 2 n+2m + {2rn - 2)c 1 fi 2n +2m-2 H h 2c m _i/i 2n +2 = (2n + l)/i 2n , (11) 

for n = 0, and /i 2j+ i = 0. /i 2 , /i4, . . . , /i 2m _ 2 are given by /i 2j //i = 



Q(x 2 )). 



4. Solution of Freud's equations. 



(J 2r ) 



The following lemma will be used: 

Lemma 1. Let = {a n ^}^ = i, k = 0, 1, ... , be the solutions of the 
quotient- difference equations 





then, F n (a,( ') satisfies 




i {F n (aW) - F n (a^)) = (F n+1 (a^) - F^a 



(k+i) 



n even 




n odd. 



(13) 



(k)\ 

or ') 




' 


0-1, k 


0-1, k 
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P'(J k )) nn _ 1 , ond where J \ is the Jacobi matrix 



0,2, k 

a 3ifc 



The quotient-difference equations are 
well known in orthogonal polynomials 
identities investigations, see [74, 75] for 
recent examples. The notation chosen 
here is motivated by the interpretation 
of a Ujk as recurrence coefficient of orthogonal polynomials related to the weight 
x 2k w(x), but we do not need this interpretation in order to establish the 
theorem and the lemma. Actually, the theorem 2 is meant to recover an 
interpretation of the solutions of (6) in terms of orthogonal polynomials, 
but arbitrary initial data a i} a 2 , . . . , a m -i will yield formal orthogonal poly- 
nomials through formal moments fij which, from (11), can still be written 



x J exp( — P(x)) dx } but where the support S is now a system of arcs 



in the complex plane, see [66, "Laplace's method"], [62,63]. 

Remark that, with a„ j0 = o, n and a 0jk = 0, all the a Ujk } s are completely 
determined by (12) as functions of a i} a 2 , . . . 

Proof of Lemma 1. The quotient-difference equations (12) can be 
written J k L k = L k J k+ i } where L k is the infinite lower triangular matrix 

One finds also J\ = L k R k , 
= R k L k , where R k is the 
transposed of L k (the famous 
LR relations). So 
L k R k L k J k+ i = LkJ^+i, 
L k R k L k J 2p 



0-1, k 






ai,k+i 






d2,k 






«2,fc+l 





«3,fc+l 










0-5, k 



J 3 k L k 



L k J 2 k P +i , for any odd power, so 



Jl P+1 L 



P'(J k )L h 



k — ^k^k^k" 

L k P'(J k+ i) } and we look at 



the element of indexes 



n. n 



1 



P '{J ' k))n,n-l(Lk)n-l,n-l + (P' (J k))n,n+l (Lk)n+l,n-l ~ 

(Lk)n,n-2(P'(Jk+l))n-2,n-l + (L k)n,n(P' \J k+l))n,n-l i 



t? I „k\ \ L k )n-l,n-l p ( k\ \ L k )n+l,n-l 

or fi n (a ) \-fi n+ i(a ) = 

On,k O n j r \^k 

(Lk)n,n-2 p / | (Lk)n,n p / 

Vl^ ) H ^n\d J, 

a n-l,k+l a n,k+l 

leading to (13). □ 
Proof of Theorem 2. From (13), with (6), (a<°) = a), and F (a (fc) ) = 
for all k, one finds by induction F n {a^) = n + k(l — ( — 1)™), so, using 
(7), fi(aW) = (J k P'(J k )) 0fi = 2k + 1, k = 0,1,... Let /i 2k = constant 
xafaf l a\ 2 • • • a\ k _ 1 and fi 2k +i = 0. From the LR relations, (J 2p )o,o = 

(L k L k+ i ■ ■ ■ L k+p R k+p ■ ■ ■ R k+ iR k ) 0j0 = Cli t k a l,k+1 ' ' ' a l,k+p-l = t l 2k+2p/ t l 2k- 
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Finally, 2k + 1 = (J k P'(J k )) 0fi = 2mc (J 2 k m ) 0fl + (2m - 2) Cl (J 2 k m - 2 ) 0fl + 

••• + 2c m _i(J^) ,o = 2mco/i 2 fc+2m//"2fc + (2m - 2)c 1 fi 2 k+2m-2/ ^2k + ••• + 
2c m _i/i 2 fc+2//W2fc, whence (11), and the determinant ratios (3) follow from 
(12): a 2 nk = H n _i^H n+ i^ I (H n ^) 2 where H Ujk is the Hankel determinant of 

rOWS /i 2 fc,/"2fc+l(= 0), . . . ,/i 2 fc+n-i; /"2fc+l, • • . ,/i 2 fc+rx, etc., Up to 

l^2k+n-l , ■ ■ ■ , l^2k+2n-2- D 

The isomonodromy approach looks at the linear differential equation sat- 
isfied by an orthogonal polynomial p n [5,16,17,37,38] and produces isomon- 
odromy identities of interest here. For instance, working (9) leads to a 2 n = 
(4a)- 1 / 2 Piv((4a)" 1 / 2 /3; -re/2, -re 2 /4) (Kitaev, [23-26,43]), where Pi V (t; A, B) 
is the solution of the fourth Painleve equation 

y = ^- + 3 -^ + 4ty 2 + 2(t 2 -A)y + ^ 

which remains 0(t~ l ) when t — > +oo [60, pp. 231-232] (the importance and 
relevance of this latter condition is not yet quite clear). See [1—3, 16—20, 22— 
26,39,43,44,46] for more on such connections. 

5. Freud's equations as discrete Painleve equations. 

As far as I understand the matter, discrete Painleve equations were first 
designed as clever discretisations of the genuine Painleve equations, so to keep 
interesting features of these equations, especially integrability. For instance, 
(9) when (3 ^ is a discretisation of the Painleve-I equation [23-26, 33, 34, 
43]. Then, these features are examined on various discrete equations not 
necessarily associated to Painleve equations. So, for instance, [33, p. 350] 
encounter (9) with (3 = as a discrete equation which is no more linked to 
the discretisation of a continuous Painleve equation (and the authors of [33] 
call (9) a discrete Painleve-0 equation when (3 = 0). 

If we want to check that the Freud's equations are a valuable instance of 
discrete Painleve equation, the following features are of interest, according to 
experts [48]: 

5.1. Analyticity. Each a 2 l+m _ 1 is a rational function of preceding elements, 
so a meromorphic function of initial conditions. 

5.2. Reversibility. There is no way to distinguish past and future, (6) is 
left unchanged when a n+ i, a n+2 , . . . are permuted with a n _i, a n _ 2 , . . . Indeed 
( [56, p.369]), (6) is a sum of terms a 2 n a 2 nJrii a 2 nJri2 ■ ■ ■ a 2 n+ip , provided ^ «i + l ^ 
Z2 + 2^ ••• ^ i p + p ^ p + I. Then, {— i i} — i 2} . . . , — i p } satisfies the same 
conditions and is therefore present too: ^ — i p + 1 ^ —i P -i + 2 ^ • • • ^ 
-i\ + p < p + 1- 

5.3. Symmetry. An even stronger property is the following: 

Theorem 3. The matrix of derivatives ( — — — - — - J is symmetric. 

\dlog a m J 
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Proof: see [56, 57]. 

5.4. Integrability. The existence of the formula (3) according to Theorem 2 
shows how (6) can be solved. The similarity of (3) with known solutions of 
discrete Painleve equations is also striking [42,45]. The structure of the 
solution leads also to 

5.5. Singularity confinement. The movable poles property of continuous 
equations is replaced by the following [34]: if initial values are such that 
some a n , say a no becomes very large, only a finite number of neighbours 
a no+ i, . . . , a no+p should be liable to be very large too; moreover, a no+p+ i } etc. 
should be continuous functions of a no _i,a no _ 2 etc. The formula (3) shows 
that the singularities of the zeros of H n which are determinants built 
with solutions of (11), therefore continuous functions of the initial data. 



6. Generalizations of Freud's equations: semiclassical 
orthogonal polynomials, etc. 

Let the weight function satisfy the differential equation w 1 /w = 2V/W } 
where V and W are polynomials. The related orthogonal polynomials on a 
support S are then called semi-classical [6-8,40,41,64,65,72] (classical: de- 
grees of V and W ^ 1 and 2). One can still build equations for the recurrence 

coefficients by working out the integrals I n = W(x)p' n (x)p n _i(x)w(x) dx, 

Js 

and J n = W(x)p' n (x) p n (x)w(x) dx. Many equivalent forms can be found, 
Js 

and it not yet clear to know which one is the most convenient. If we manage to 
have W(x)w(x) = at the endpoints of the support S, explicit formulas fol- 
low [77] (see also formula (4.5) of [1]). For instance, with w(x) = exp(— Ax 2 ) 
on S = [—a, a] (see [15] for other equations), one takes W^z) = x 2 — a 2 
and V(x) = —Ax(x 2 — a 2 ), uses p' n = np n _\ja n + [2{a\ + • • • + a 2 l _ 1 ) — 



na 



2 1 ]p n _ 3 /(a n _ 2 a n _ia n ) + • • • from (2), to find 2(a\ + • • • + a 2 ^) + (2n + 



l)a 2 n - no 2 - 2Aa 2 n (a 2 n _ 1 + a 2 n + a 2 n+1 - a) = 0, or (2n + l)x n+1 - (2n - l)x n - 
na 2 - 2A(x n+ i - x n )(x n+2 - x n _ x - a) = 0, where x n = a\ + • • • + a 2 n _ x . 

It is not clear how to recover reversibility and symmetry. For an even 
much more nasty case, see [61]. 

But the root of the matter lies even deeper, as already recognized by La- 
guerre in 1885 [49]! The point is to derive equations for the coefficients in 

the Jacobi continued fraction f(z) = from a differen- 

a{ 

z - b 

z — b\ — ■ ■ ■ 

tial equation Wf = 2V f + U with polynomial coefficients (same W and 
V as above: orthogonal polynomials lead to this continued fraction through 

f(z) = (z — x)~ l w(x) dx). See [30,72] for this technique. 
Js 
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The continued fraction approach leads readily to the Riccati extension 
Wf = Rf 2 + 2Vf + U [20,55], as each continued fraction 

f k (z) = 2 1S f° un( i to satisfy a Riccati equation Wf' k = 



z - b k+1 

— a^0/ ; _i/|+2O/ ; // ; — 0/; too, where 0^ and ft k are polynomials of fixed degrees 
(i.e., independent, of k). Equations for the recurrence coefficients a k s and 
b k s follow from recurrence relations for the ©fc's and fi^'s (use fo(z) = f(z) 
and f k (z) = l/[z - b k - al +1 f k+1 (z)]). 

Finally, the differential operator can be replaced by a suitable difference 

operator, amounting to W{z) ~ = R ( z ) f M z )) f (^{ z )) 

-\-V(z)[f((fi(z)) + f((f 2 (z))] + U(z), where fi(z) and (f2(z) are the two roots 
of a quadratic equation A(f 2 (z) + 2Bz<p(z) + Cz 2 + 2D(f(z) + 2Ez + F = 0, 
see [36,58]. 
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26: In felaweshipe, and pilgrimes 
were they alle, 

27: That toward caunterbury 
wolden ryde. 

28: The chambres and the stables 
weren wyde, 

29: And wel we weren esed atte 
beste. 



In fellowship, and pilgrims were 
they all 

That toward Canterbury town 
would ride. 

The rooms and stables spacious 
were and wide, 

And well we there were eased, and 
of the best. 
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